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Using 467 fb~^ of e"'"e~ annihilation data collected with the BABAR detector, we measure 
B(r--^^^,-y'^.) ^(0.9796 ± 0.0016 ± 0.0036), SC^'^^^-^^) =(0.5945 ± 0.0014 ± 0.0061), and 

B(T-^^Vell) =(0-03882 ± 0.00032 ± 0.00057), where the uncertainties are statistical and system- 
atic, respectively. Prom these precision r measurements, we test the Standard Model assumption of 
li-e and t-h charge current lepton universality and provide determinations of \Vua\ experimentally 
independent of the decay of a kaon. 



PACS numbers; 11.30.Hv, 12.15.Hh, 13.35.Dx, 14.60.Fg, 14.40.Aq, 13.66.Lm 
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Decays of the r lepton to a single charged particle and 
neutrino(s) probe the Standard Model (SM) predictions 
of charged current lepton universality and the unitar- 
ity relation of the first row of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix jlj. Previous mea- 
surements of universality expressible in terms of the 
coupling strength (g^) of lepton of flavor £ to the charged 
gauge boson of the electroweak interaction are in agree- 
ment with the SM where Qr/Pfi = — 1- Similarly, 
kaon decay measurements [3|,|j] sensitive to \ Vus\, the rel- 
ative weak coupling between up and strange quarks, yield 
a value consistent with unitarity ( | Vud P + 1 ^usP + 1 Vub \ ^ = 
1) where nuclear beta decays provide \Vud\ M \Vub\ 
IS negng ible faj. However, new physics that couples pri- 
marily to the third generation could be revealed through 
deviations from the SM in precision universality and |14s | 
measurements involving the r. Significant deviations of 
this nature are unambiguous signatures of new physics 
that provide crucial but complimentary information to 
the direct searches for Higgs bosons [6|and other new 
physics models with e.g. lepto-quarks [7|, heavy gauge 
W' or Z' bosons, heavy quarks or leptons, composite- 
ness or extra dimensions [8|. 

Recent measurements of the sum of strange r branch- 
ing fractions interpreted in the framework of the Oper- 
ator Product Expansion (OPE) and finite energy sum 
rules yield a value of \Vus\ that is approximately three 
standard deviations (cr) lower than expectations from 
CKM unitarity This paper addresses both ex- 

perimental and theoretical aspects of this question by 
providing the first precision measurements of Rk = 



B(r-^e-i7,.,) 110| and Rk/. = B{r-Z.~Z) enabled by 
the unique combination of a very large r sample with par- 
ticle momenta amenable to particle identification using 
Cherenkov radiation. By using values of the meson decay 
constants from lattice QCD , we provide two precision 
determinations of \Vus \ from r decays independent of the 
OPE framework. We also report on new measurements 
or Ut, ana it^ _ ., ^ . k„ 



vides an improved measurement of t/^/ge whereas i?^ and 
Rki when compared to the muonic branching fractions 
of the pion and kaon, yield improved measurements of 
dr/gti involving pseudoscalar mesons. 

The data sample corresponds to an integrated lumi- 
nosity of £ = 467 fb^^ recorded at an e'^e~ center-of- 
mass (CM) energy (y/s) near 10.58 GeV and was col- 
lected with the BaBAR detector at the SLAC PEP-11 
e+e" storage rings. With a luminosity- weighted av- 
erage cross-section of CTe+e-^T + T- = (0.919 ± 0.003) 
nb 13, l3|) this corresponds to the production of 4.29 x 



10^ T-pair events. The BaBAR detector [13j is com- 
posed of a silicon vertex tracker, drift chamber (DCH), 
ring-imaging Cherenkov detector (DIRC), and electro- 
magnetic calorimeter (EMC), all contained in a 1.5-T 
solenoid. The iron flux return for the solenoid is instru- 



mented (IFR) to identify muons. 

Tau-pair events are simulated with the KK Monte Carlo 
(MC) generator JJj], which includes higher-order radia- 
tive corrections. We simulate r decays with TAUDLA [15| 
and PHOTOS [l6[ using measured branching fractionsjs]. 
The detector response is simulated with GEANT4 [l!l]. 
Simulated events for signal as well as background pro- 
13- is!] are reconstructed in the same manner as 



cesses 
data. 



The MC samples are used for selection optimiza- 
tion, control sample studies, and systematic error studies. 
The number of simulated non-signal events is comparable 
to the number expected in the data, with the exception of 
Bhabha and two-photon events, which are not simulated 
but which data studies show to be negligible. 

We study e+e^ — >■ r+r^ events with the t~ decaying 

via — > e^VeVr, pi- Vf^Vr, — > 11^ Vt Or — > 

K^Vr modes and the decaying via a — >■ tt+tt+tt"!/,- 
tagging channel with the selection criteria optimized to 
minimize the combined statistical and systematic un- 
certainties jl^l- The number of signal events for decay 
modes i = {e, [i, tt, if} = {e~Vf,VT,\i~V ^v^^ t^~Vt-, K~i'r} 
are Nf = (^N^ — iVf ) where £i is the efficiency (in- 
cluding B{t- ^ TT-TT-TT+Ur) = (8.85 ± 0.13)% Q ) , TVf 

the number of selected data events, and the estimated 
number of background events for the i*^ mode. 

We measure the ratios Ri = Nf/Nf which normalizes 
to the most precisely known relevant SM process avail- 
able, and in which several common sources of systematic 
uncertainity cancel. are multiplied with reproducible 
random numbers until all efficency and uncertainity es- 
timates are finalized. Once unblinded, we use the values 
of the three branching ratios to update world averages of 
the branching fractions, which we then use to recalculate 
the backgrounds for our final results. 

Events with a net charge of zero and with four well- 
reconstructed tracks not originating from the conversion 
of a photon in the detector material are selected. For 
good particle identification, each track is required to be 
within the acceptance of the DIRC and EMC, and have 
a transverse momentum greater than 0.25 GeV to ensure 
that it reaches the DIRC. The plane normal to the thrust 
axis divides the event into hemispheres in the CM frame. 
The "signal" hemisphere contains a single track and the 
"tag" hemisphere the other three tracks. 

Each tag hemisphere track is required to be consis- 
tent with being a pion and the energy deposited in the 
EMC unassociated with any tracks in this hemisphere is 
required to be less than 0.20 GeV. Also, events that con- 
tain track pairs consistent with coming from a Kg are 
vetoed. 

The signal track momentum is required to lie between 
1 and 4GeV/c. Information from the five detector sub- 
systems is combined in likelihood selectors which iden- 
tify e, TT, and K particles and in a neural network which 
identifies muons. The n-K separation is provided by the 
DIRC and DCH whereas 7r-/i separation is primarily ac- 
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complished with the IFR and EMC. The identification 
efRciencies are given in Table U and cross-contaminations 
are given below. We suppress di-muon and Bhabha back- 
grounds by requiring signal tracks identified as a lepton to 
have CM momentum less than 80% of ^/sjlc. To reduce 
cross-feed from e into the tt and K channels, the ratio of 
deposited electromagnetic energy of a tt or K candidate 
track to its measured momentum, E/pc, is required to 
be less than 0.85. A pion track also passing a loose muon 
selection is rejected. A similar veto is applied for a kaon 
track passing the loose muon selection if its measured 
momentum exceeds 3 GeV/c. Also, events with an EMC 
energy > {1.0, 0.5,0.2,0.2} GeV in the signal hemisphere 
unassociated with the {e, /i, tt, K} track are removed. 

Pion and kaon control samples from — > 
7r+Z)°,D° — > t:^K~ decays are used to study and cor- 
rect for small differences between MC and data. We 
cross-check these with independent 7r+ [K^] control 
samples from t~ — >■ ■n~ 'k~ tt'^ Vr (r^ K^'k^K^Vt) 
decays using particle identification of two of the oppo- 
sitely charged particles and the fact that the wrong sign 
— > Ti^ TT~ K'^Vt decays are heavily suppressed. Sam- 
ples of radiative Bhabha and radiative /x-pair events pro- 
vide control samples of electrons and muons. The system- 
atic uncertainty associated with charged particle iden- 
tification is assessed from the control sample statistical 
errors, consistency between control samples, and the sen- 
sitivity of the control sample corrections to the number 
of particles near the track. The statistical errors in the 
more limited cross-check control samples dominate these 
errors. Because we use control samples to correct charge 
conjugate particles separately, charge-dependent detec- 
tor responses are accounted for by construction. 

To remove two-photon and Bhabha backgrounds, the 
event must have a missing CM energy between 10% and 
70% of -yi. The angle between the missing momentum 
and electron beam direction in the CM, 6*^^^, is con- 
strained to satisfy |cos(6'^?fJ| < 0.7, the thrust of the 
event is required to be above 0.9, and the net missing 
transverse momentum in the CM greater than 0.009-y/s/c. 

Each of the three tag-side tracks has an electron veto 
applied to further reduce the Bhabha contamination. 
This results in less than 0.03% contamination from two- 
photon events and less than 0.1% contamination from 
Bhabha events in the electron signal sample. These back- 
grounds were investigated by studying samples enriched 
in Bhabha and two-photon events by adjusting the re- 
quirements on the thrust, cos(6'^5fg), and transverse mo- 
mentum of the event. Potential background from Bhabha 
events were further probed by studying the number of 
events having a high signal track momentum as the elec- 
tron veto was progressively lifted from one, then two, and 
finally all three tracks in the tag hemisphere. 

To suppress backgrounds in the — > ti^Vt and 
— > K~Vr channels from r decays with undetected 
neutral particles other than the Vr (e.g. K'^ mesons, v^), 
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FIG. 1: Data (points) and MC (histograms) distributions of 
CM momentum for (a) t~ e^VeVr, (b) r~ fi~Vf^.Vr, 
(c) T~ — >■ 7r~z/T- and (d) — > K~Vr modes. The small 
differences between MC and data are accounted for in the 
systematic errors. 



we reconstruct the direction of the back-to-back t+t~ 
system in the CM frame. The polar angle of the r mo- 
mentum with respect to the tag-side hadronic system is 
calculated assuming that the CM energy of the r is -\/s/2, 
and the azimuthal angle of the r momentum is fixed to a 
value that has been optimized to minimize the total error 
on Bk/tt (20j . With this estimator for the r momentum, 
we require the missing mass in the signal hemisphere to 
be less than 0.56 GeV/c^. 

For the selected t~ fi~T'f^i'r events, the domi- 
nant backgrounds are — > tt^Vt (1.46 ± 0.01)% and 
T- TT-TT°Ur (0.85 ± 0.01)%. For the tt^i/^ 
channel, the dominant backgrounds are — > ^^V^j^Vr 
(12.90 ±0.07)%, T- ^ TT-TT^Vr (5.87 ± 0.04)%, and non- 
T backgrounds (0.34 ± 0.05)%. The major backgrounds 
in the t" K~'Vr channel are from r" — >■ ■K~'Vr decays 
(10.06 ±0.13)%, ^ K-Klvr (3.87 ±0.41)%, t" ^ 
K-ii°Vr (1.97 ± 0.14)%, T- -K-TT^Vr (1.07 ± 0.06)%, 
and non-r backgrounds (2.58 ± 0.38)%. The uncertain- 
ties are from MC statistics, branching fractions and, for 
non-T backgrounds, the systematic uncertainty on back- 
ground rates. Fig. [T] shows the momentum distributions 
in the CM frame for each of the four decay modes for 
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data, along with the background MC contributions. 

For the t~ — >■ e~Ve.v.T channel, 884426 events are se- 
lected with an efficiency and purity of (0.589 ± 0.010)% 
and (99.69±0.06)%, respectively. The number of selected 
events, efficiency, purity and systematic uncertainties on 
Ri of the T" ^~T>^Vt, t~ — ?► t^'Vt, and — > K^Vr 
selections are presented in Table HI These uncertain- 
ties include contributions from the particle identification, 
the sensitivity to detector response including the impact 
of changing the MC momentum scale and DCH resolu- 
tion, modelling of hadronic and electromagnetic showers 
in the EMC, the EMC energy scale, and angular mea- 
surements made by these detectors within their mod- 
elling uncertainties, the backgrounds, initial- and final- 
state radiation, radiation in r decays, rate and shape of 
T~ 7r~7r~7r+^'T- decays, the trigger, and Ca^+f.-^r+T- ■ 
The systematic uncertainty on is dominated by uncer- 
tainties in particle identification. The Rt^ and Rk mea- 
surements have additional dominant contributions from 
the detector modelling and associated backgrounds, due 
to stronger cuts on the EMC energy necessary to reduce 
non-T backgrounds. Presence of the ~20% backgrounds 
in these channels render them more sensitive to the mod- 
elling of the tag-side decays. The dominant background 
uncertainty in the Rt^ measurement arises from the elec- 
tron contamination in the tt sample investigated by mea- 
suring the number of events that fail the E/p electron 
veto requirement in data and MC. In the Rk event sam- 
ple, the uncertainty arising from r decay branching frac- 
tions of background modes is 0.58%, which is dominated 
by the uncertainty of the r~ — > K^K~Vr fraction. There 
is also a 0.49% uncertainty assigned for qq backgrounds, 
which are studied using events with an invariant mass of 
the tracks in the tag hemisphere above the T-mass and 
cross-checked in regions of thrust and cos(6'^^j,) enriched 
with these backgrounds. 



TABLE I: Number of selected events, purity, total efficiency, 
component of the efficiency from particle identification, and 
systematic uncertainties (in %) on Ri for each decay mode. 







TT 


K 




731102 


369091 


25123 


Purity 


97.3% 


78.7% 


76.6% 


Total Efficiency 


0.485% 


0.324% 


0.330% 


Particle ID Efficiency 


74.5% 


74.6% 


84.6% 


Systematic uncertainties: 






Particle ID 


0.32 


0.51 


0.94 


Detector response 


0.08 


0.64 


0.54 


Backgrounds 


0.08 


0.44 


0.85 


Trigger 


0.10 


0.10 


0.10 


■K~-K~-K^ modelling 


0.01 


0.07 


0.27 


Radiation 


0.04 


0.10 


0.04 


B{t~ — ^ -k' ■K~ -K^ Ut) 


0.05 


0.15 


0.40 




0.02 


0.39 


0.20 


Total [%] 


0.36 


1.0 


1.5 



The measured branching ratios and fractions are: 
R 



Rti 

Rk 

Rh = Rtt + Rk 
Bir- K-Vr) 



(0.9796 ±0.0016 ±0.0036) 
(0.5945 ±0.0014 ±0.0061) 
(0.03882 ± 0.00032 ± 0.00057) 
(0.6333 ±0.0014 ± 0.0061) 
(17.46 ±0.03 ±0.08)% 
(10.59 ±0.03 ±0.11)% 
(0.692 ±0.006 ±0.010)% 



where h — t: or K and we use B{t — e Vei^r) = 
(17.82 ± 0.05)% [3]. The off-diagonal elements of the 
correlation matrix for the measured ratios (branching 
fractions) are =0.25 (0.34), =0.12 (0.20), and 
PttK =0.33 (0.36). The fi and tt measurements are con- 
sistent with and of comparable precision as the world av- 
erages Q whereas the K measurement is consistent with 
but twice as precise as the world average fsl. 
Tests of /X — e universality can be expressed as 



B{t fi v^Vr) f (ml/ml) 



B{t- 



VeVr) f (ml/ml)' 



where f{x) = 1 — 8a; + 8a;'^ — a;"* — 12a;^logx, assuming 
that the neutrino masses are negligible |2l| . This gives 

i^Y"^ ~ 1.0036 ± 0.0020, yielding a new world average 
of 1.0018 ± 0.0014, which is consistent with the SM and 
the value of 1.0021 ± 0.0015 from pion decays [3, [III. 
Tau-muon universality is tested with 



B{t hi^r) '2mhmf^Th / 1 - mf^/ml 
B{h yLV^_,) (1 + Sh)m^Tr I 1 - m\/ml 



where the radiative corrections are 5.^ = (0.16 ± 0.14)% 
and 5k = (0.90 ± 0.22)% [23]. Using the world averaged 
mass and lifetime values and meson decay rates Q, we 

determine = 0.9856 ± 0.0057 (0.9827 ± 0.0086) 

and ( — I = 0.9850 ± 0.0054 when combining these re- 

suits; this is 2.8(7 below the SM expectation and within 
2ct of the world average. 

We use the kaon decay constant Jk — 157±2 MeV |11| , 
and our value of 



us\ Tfl T-j- 



1 



S 



EWi 



where Sew = 1.0201 ± 0.0003 [2j|, to determine = 
0.2193 ± 0.0032. This measurement is within 2cr of the 
value of 0.2255 ±0.0010 predicted by CKM unitarity and 
is also consistent with the value of = 0.2165±0.0027 
derived from the inclusive sum of strange r decays [§] . 

Both of our measured \yus\ values depend on absolute 
strange decay rates. Our value of Rki-k = (0.06531 ± 
0.00056±0. 00093), however, provides a value driven 
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by the ratio between strange and non-strange decays. We 
use Jk/U = 1-189 ± 0.007 [111, |Kd| i, and the long- 
distance correction S^o = (0.03 ± 0.44)% estimated [25| 
using corrections to t — > and h ^ iiv [2^ [2^ in 



Pk\Vu. 



(I + Sld), 



to obtain |Ks| = 0.2255 ± 0.0024 where short-distance 
electro-weak corrections cancel in this ratio. This value 
is consistent with CKM unitarity [5] and 2.5(7 higher than 
\Vus\ from the inclusive sum of strange r decays. 
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